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ABSTRACT 

We report on the discovery of the z = 1.016 cluster RzCS 052 using a modified 
red sequence method, foUowup spectroscopy and X-ray imaging. This cluster has a 
velocity dispersion of 710 ± 150 km s"'^, a virial mass of 4.0 x 10^^ Mq (based on 
21 spectroscopically confirmed members) and an X-ray luminosity of (0.68 ± 0.47) x 
ID*** ergs s~^ in the [1-4] keV band. This optically selected cluster appears to be 
of richness class 3 and to follow the known Lx — cry relation for high redshift X- 
ray selected clusters. Using these data, we find that the halo occupation number for 
this cluster is only marginally consistent with what expected assuming a self-similar 
evolution of cluster scaling relations, suggesting perhaps a break of them at 2 1. 
We also rule out a strong galaxy merging activity between z = 1 and today. Finally, 
we present a Bayesian approach to measuring cluster velocity dispersions and X-ray 
luminosities in the presence of a background: we critically reanalyze recent claims for 
X-ray underluminous clusters using these techniques and find that the clusters can be 
accommodated within the existing Lx — cry relation. 

Key words: Galaxies: evolution — galaxies: clusters: general — galaxies: clusters: 
individual RzCS 052, — (Cosmology:) dark matter — X-rays: galaxies: clusters — 
Methods: statistical — Galaxies: kinematics and dynamics 



1 INTRODUCTION 

Clusters of galaxies are not only a powerful tool to study 
galaxy evolution but can also be used to constrain cosmo- 
logical parameters, resolving several parameter degeneracies 
(e.g., Allen et al. 2004; Albrecht et al. 2006). In particular, 
clusters at high redshifts {z > 1), of which only a handful 
are currently known, provide the greatest leverage in deter- 
mining the nature of the acceleration constant (e.g., Rapetti 
2007). These determinations, however, rely on an accurate 
estimate of the cluster mass, whose uncertainty is arguably 
the dominant contributor to the error budget in deriving 
cosmological parameters from cluster statistics (Henry 2004; 
Albrecht et al. 2006). 

Ideally, one wish to apply the virial theorem to get a di- 
rect measurement of cluster masses. In fact, the dark matter 
velocity dispersion is an extremely good tracer of the halo 
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masses in all simulations (Evrard et al. 2007), and galax- 
ies are nearly unbiased velocity tracer (Evrard et al. 2007 
and references therein; Rines Diaferio & Natarajan 2007), 
in good agreement previous works (Biviano et al. 2006, Tor- 
men et al. 1997). The measurement of the cluster velocity 
dispersion requires a large number of radial velocities, which 
are observationally expensive to obtain, especially for high 
redshift clusters. For this reason and because each mass es- 
timator carries some key informations, more commonly the 
scaling between pairs of more easily observable mass-related 
quantities is studied, such as X-ray luminosity, temperature 
or the Yx (Kravtsov et al. 2006) parameter, or optical rich- 
ness. These studies often look for outliers, however their 
search is blessed by data limitation: for example in the search 
of clusters X-ray dim for their optical richness, Donahue et 
al. (2001) and Gilbank et al. (2004), both mostly worked 
with putative clusters (i.e. not spectroscopically confirmed) 
and X-ray undetections. 

Only few works directly address the relative quality of 
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Figure 1. True-color (z' [3.6] [4.5]) degraded-resolution (to make galaxies not too small when printed) image of a region of a 24 Mpc'^ 
area around RzCS 052. Spectroscopically confirmed clusters and isodensity contours for red galaxies are also marked. Note the number 
density contrast of reddish galaxies between the cluster center and the right part of the image. The ruler is 1 arcmin long; North is up, 
East is to the left. 



different mass estimators with velocity dispersion: Borgani 
& Guzzo (2001) compare the scatter of two mass estimators, 
X-ray luminosity and richness, and found that the former is 
a better mass tracer than the latter when the former is uni- 
formly measured and the latter is taken from a 50 years old 
paper reporting eye-estimate of the cluster optical richness 
(the Abeh 1958 catalog). In both CNOC and nearby clus- 
ters, mass correlates better with richness than with X-ray 
luminosity (Yee & Elhnson 2003; Popesso et al. 2005). Eke 
et al. (2004) found that optical luminosity is a better proxy 
of mass than velocity dispersion in common conditions, i.e. 
when velocities are available for a small sample of galaxies. 
A related issue, which we will examine below, is whether 
it exists clusters X-ray dim for their mass (velocity disper- 
sion), (e.g. Lubin, Mulchaey & Postman 2004, Fang et al. 
2007, Johnson et al. 2006). 

The relation between richness and mass has received 
some recent attention in the form of the halo occupation 
function (Berlind & Weinberg 2002; Lin et al. 2004 and ref- 
erences therein) whose first moment is the halo occupation 
number (HON) , the average number A*' of galaxies per clus- 
ter of mass M. In order to address the evolution of the 
HON, velocity dispersion information is often unavailable 
for a large cluster sample, mass and cluster size are inferred 
from other mass-related quantities (for example the X-ray 
temperature), and assumed to evolve self similarly. The evo- 
lution of the HON with redshift is still unclear: the initial 
study by Lin et al. (2004) claimed that the HON increases 
at high redshift, but Lin et al. (2006) find evidence that it 
does not evolve strongly out to z ~ 1, suggesting that the 
galaxy population in clusters was established and assembled 
at early epochs. Muzzin et al. (2007) confirms the above, 
with a sample of reduced redshift leverage and hence re- 



duced evolution sensitivity, but available velocity dispersion 
information. 

Here, we present the photometric discovery, spectro- 
scopic confirmation and X-ray properties of a new z — 1.016 
cluster of galaxies (RzCS 052), a cluster optically rich but 
undetected in the XMM-LSS survey (Pierre et al. 2007) , and 
hence possibly X-ray dark (i.e. dim for its mass). We derive 
its global properties (richness. X-ray luminosity, velocity dis- 
persion and mass) and study these in the context of cluster 
scaling relations {Lx — cr, HON) at high redshift. In par- 
ticular, we test the claim that the HON (the way galaxies 
populate cluster-scale haloes) has not changed z ~ 1 (Lin et 
al. 2006) under far less assumptions than the original claim. 
We also present a Bayesian approach to the determination 
of cluster velocity dispersion and X-ray luminosity and use 
it to critically examine recent claims about the existence of 
underluminous X-ray clusters. A companion paper (Andreon 
et al. 2007) addresses the use of RzCS 052 as a laboratory 
for studying galaxy formation and evolution. 

We adopt = 0.7, = 0.3 and Ho = 70 km s"^ 
Mpc~^. Magnitudes are quoted in their native photometric 
system (Vega for R, SDSS for z'). 



2 THE DATA & DATA REDUCTION 

2.1 Photometry: CTIO Rz' images 

Broadband images for a 7 deg^ region around this cluster 
were obtained at the Cerro Tololo Inter-American Obser- 
vatory Blanco 4m telescope in the R and z' (Ac ~ 9000A) 
filters using the Mosaic II camera. We use the same imaging 
data as Andreon et al. (2004a), where details on the data and 
their analysis are found. Briefiy, Mosaic II is a 8192 x 8192 
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Figure 2. Spectra of RzCS 052 members coming the VLT run. We have vertically shifted the spectra and zoomed on a reduced wavelength 
range for display purpose. 
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Figure 3. Left: Colour-magnitude diagram for galaxies (close circles) within 1 arcmin from the cluster center or with a known redshift 
(red closed circles for members, blue crosses for interlopers) using CTIO discovery data. R and z' mag completeness limits are show 
with dashed lines. The green line is the expected CM at the cluster redshift, from Kodama & Arimoto (1997) Right: Colour histograms 
of galaxies within 2 arcmin from the cluster center (solid histogram), and of the average control field (measured on a 0.36 deg^ area, 
dashed histogram), normalized to the cluster area. A clear excess is seen, especially at ij — 2 > 1.5 mag. The hashed histogram is the 
colour distribution of spectroscopically confirmed member galaxies. In both panels, a few objects with spectroscopic redshift are missing 
because their fall on bad CCD regions or have extreme colours. 



CCD camera with a 36' field of view at prime focus. Ex- 
posure times were 1200s in R and 1500s in z': seeing was 
between 0.9" and 1.0" in the final images. 

Source detection and photometry were carried out us- 
ing SExtractor (Bertin & Arnouts 1996). Colours and mag- 
nitudes were computed within a fixed 2" (radius) aperture, 
and corrected for minor differences in seeing, as in Andreon 
et al. (2004a). Completeness magnitudes {5a in a 3" aper- 
ture), computed as in Garilli, Maccagni & Andreon (1999) 
are R = 24.0, z' = 22.5 mag. 

Figure 1 shows a true-colour image of RzCS 052, as 
derived from CTIO z' image and IRAC Spitzer [3.6] and 



[4.5] images. Spitzer data reduction is described in Andreon 
(2006a), that also presents the composite stellar mass func- 
tion and the 3.6 /im luminosity function of many clusters, 
including RzCS 052. 



2.2 Spectroscopy 

Multiobject spectroscopy was carried out on Gemini in late 
2003, and VLT in late 2003 and during 2004. On VLT, the 
spectra were taken using FORS2 with the GRISJ500I and 
the OG590 filter for a total integration time of 11 ks. On 
Gemini, the spectra were obtained with the Gemini Multi 
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Table 1. J2000 coordinates and redshift of galaxies within 4000 
km s^i (rest-frame) of RzCS 052 
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Object Spectrograph (CMOS), operating in nod & shuffle 
mode (Abraham et al 2004; Cuillandre et aL 1994) in order 
to perform accurate sky subtraction, with the R150 grating 
for a total integration time of 15 ks. 

The Gemini GMOS package for IRAF was used to cal- 
culate the wavelength solutions and to reduce the multiob- 
ject observations into one-dimensional spectra. The RVSAO 
package (Kurtz & Mink 1998) was used to measure redshifts 
(and their errors) of target galaxies by cross-correlation with 
stellar and galaxy templates of known radial velocity (Tonry 
& Davis 1979). 

A total of 57 spectra of 54 galaxies yielded reliable red- 
shifts, with typical individual errors on redshift of 50 to 150 
km/s (depending on instrument, exposure time, object spec- 
trum, etc.). Three galaxies with duplicate observations have 
concordant redshifts in the two data sets. Figure 2 shows 
the spectra of RzCS 052 members from the VLT run. Table 
1 list position and redshift of galaxies within 4000 km 
of RzCS 052. 



2.3 XMM-Epic data 

RzCS 052 was observed with XMM Newton using the Euro- 
pean Photon Imaging Camera (EPIC) instrument (Jansen 
et al. 2001) in 2002 in full-frame mode with the thin filter. 
After flares filtering, the good exposure time is ~13 ks for 
the MOSes (Turner et al. 2001), and a ~8 ks for the PN 
(Struder et al. 2001). By using the XMM-Newton Science 
Analysis System (SAS, v. 7) package and our own scripts, 
we kept only patterns between and 12 for MOS and to 4 
for PN. We flagged bad pixels, bad columns and CCD gaps, 
regions not seen by all three instruments, as well pixels con- 
taminated by the flux of other sources. We remove the en- 
ergy band [0.60-0.70] keV, where an instrumental line shows 
up because this flattens the sky background, and hence de- 
creases the complexity of the model used to describe its 
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Figure 4. Redshift distribution of all successfully measured red- 
shift in the cluster line of sight (top panel) and around the 
z = 1.016 (bottom panel). In the bottom panel, the curve mark 
the mean model, and the gray (yellow, in colour) region is the 68 
% highest density posterior interval. 



spatial distribution. We merged the three instruments to 
improve S/N. 

For comparison, we also reduced EPIC observations of 
a cluster at almost identical redshift, XLSSC 029 nt z — 1.05 
(Andreon et al. 2005), just one degree apart from RzCS 052. 
To make the comparison straightforward we cut the XLSSC 
029 exposure to match (almost) exactly the exposure time 
of RzCS 052. 



3 RESULTS 

3.1 Photometric discovery and colour— magnitudes 
relation 

RzCS 052 was been initially detected in 2000 using pho- 
tometric data (Rz') as a clustering of sources of similar 
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Figure 5. Spatial distribution of red galaxies (open points), spec- 
troscopically confirmed members (solid red points), and spectro- 
scopically confirmed inderlopers (blue crosses). 
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Figure 7. Radial profile of the foreground X-ray point source 
(black solid histogram) and of the XMM PN PSF (red continuous 
line), at the off-axis angle of the source, and a /3 = 2/3 model 
having rc = 16 arcsec (150 kpc if at the cluster distance, dashed 
blue line). The source is unresolved at the XMM resolution, the 
slight excess with respect to the PSF being due to the unremoved 
contribution of RzCS052, and far more compact than a /3 model 
with a typical cluster core radius at 2: = 1 . 



colour using our own version (Andreon 2003; Andreon et 
al. 2004a,b) of the red-sequence method (Gladders & Yee 
2000). This is shown in the right-hand side panel of Fig- 
ure 2, which shows that in a 2' circle centered on RzCS 
052 (02:21:42 -03:21:47 J2000) there are 13 galaxies redder 
than R—z' = 1.55 mag (solid histogram), while the expected 
number in the same area (i.e. background, average number 
measured in a 0.36 deg'^ area all around the cluster) is ~ 2.1, 
a very obvious overdensity detection. 

The colour of the red sequence (Figure 2) , compared to 
those of other high redshift clusters presented in Andreon 
et al. (2004a) suggests a redshift of « ~ 1.0. F. Barrientos 
(2006, private communication) confirmed that this cluster 
has also been detected by their red-sequence cluster survey 
with which we share the CTIO imaging. 

The right panel of Fig 2 shows that the colour distribu- 
tion of RzCS 052 galaxies (the area between the solid and 
dashed histogram) is bimodal, displaying a narrow peak at 
R~ z' — 1.6 mag and a broad excess at bluer colours. From 
now on, we define galaxies as red if 1.4 < R — z' < 1.9 mag. 
Isodensity contours for red galaxies are shown in Fig. 1 



3.2 Spectroscopical conflrmation and velocity 
dispersion 

The upper panel of Figure 4 shows the distribution of suc- 
cessfully measured redshifts in the cluster line of sight. The 
clear peak at z ~ 1.02 is in good agreement with the pho- 
tometric redshift inferred from the colour of the red se- 
quence {z ~ 1.0). The lower panel shows a detailed view 
around the cluster redshift. We measure Zduster = 1.016 
and av = 710 ± 150 km s~^ (see Appendix B for methods). 
The gapper or biweight estimators (Beers et al. 1990) give 
identical Ov 



3.3 Richness 

Figure 5 shows the spatial distribution of red galaxies (open 
points) in a wide area of 133 Mpc^ around RzCS 052. Two 
galaxy overdensities are quite obvious, both colour-detected 
by our cluster detection algorithm. Spectroscopically con- 
firmed RzCS 052 members (solid red points) and spectro- 
scopically confirmed interlopers (crosses) are also marked. 
The two large circles have a radius of 3.6', which is 80 % of 
the Abell (1958) radius (at the RzCS 052 redshift). 

We fit a /3 profile to the distribution of galaxies (Ap- 
pendix A) and remove contamination using counts from 
the field in Fig. 5, discarding the region around RzCS 058. 
Within one Abell radius we find 56 ± 20 red cluster galaxies 
brighter than z' = 22.5 mag. This number must be corrected 
to -\- 2 using the luminosity function and the 30% blue 
fraction measured in Andreon et al. (2007). The total num- 
ber of galaxies is ~ 150, that qualifies RzCS 052 as an Abell 
richness class 3. A different richness estimate is presented in 
sec 4. 

3.4 X-ray luminosity 

The RzCS 052 cluster is within the XMM Large Scale Struc- 
ture Survey (LSS) field, but not X-ray detected by the cur- 
rent XMM-LSS pipeline (Pierre et al. 2007), even thought 
several other z ^ 1 clusters are (Valtchanov et al. 2004, An- 
dreon et al. 2005, Bremer et al. 2006 and some more yet 
unpublished) . 

The left panel of Fig 6 shows the X-ray image of RzCS 
052. The X-ray source close to the optical cluster centre 
is not extended (Figure 7) and appears to coincide with a 
foreground spiral (as classified from Hubble Space Telescope 
images presented in Andreon et al. 2007) galaxy, and thus 
not associated with RzCS 052. Therefore, the flux from this 
source is discarded in the determination of the X-ray flux of 
RzCS 052. 
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Figure 6. [0.5-2.0] keV image of RzCS 052 (left panel) and XLSSC 029 (right panel), at very similar redshift and with matched X-ray 
images and smoothing. Pixels affected by other sources, or falling on CCD gaps are marked by regions. Simple eye inspection confirms 
that RzCS 052 is much fainter than XLSS 029. 



The right-hand panel of Fig. 6 shows an X-ray image of 
a very similar (in redshift, off-axis angle and exposure time 
of observations) cluster: XLSSC 029. XLSSC 029 is much 
brighter than RzCS 052. 

In order to determine the X-ray luminosities of XLSSC 
029 and RzCS052, we use the method described in Ap- 
pendix A. We assume uniform priors, zeroed in the unphysi- 
cal ranges (negative core radius, negative background inten- 
sity, negative central cluster intensity) and in ranges that 
make the total flux infinity (i.e. for (3 < 0.4). Besides re- 
turning a flux uncertainty that account for the covariance 
of all parameters, we also account for the cluster flux in the 
background region (for XLSSC 029 this turns out to lead to 
an underestimate of its count rate by 30%). 

For XLSSC 029 we found: Lx[l - 4] keV band: 4.4 ± 
0.8 10** erg s~^, formally for a temperature of 4 keV (taken 
from Pierre et al. 2006), but actually for a range of temper- 
atures because of our choice of quoting luminosities in the 
[1-4] keV band, i.e. in a band that, at the cluster redshift, 
matches the observer frame [0.5-2] band, and because of the 
very tiny dependency of the conversion factor on tempera- 
ture. For RzCS 052 we found: fx = 1.2 ± 0.8 10"^^ erg s'^ 
cm"^ and [1-4] keV rest-frame band Lx = 0.68 ± 0.47 lO** 
erg s~^ (both values are posterior mean and standard devi- 
ation). The cluster is not an ^ 1.5 detection, however: the 
posterior probability p{fx < foldata) goes to zero at fluxes 
/o > 1 lO"'^^ erg s~^ cm~^, i.e. the source has so many de- 
tected photons that data cannot be described by a model 
with a cluster signal fainter than /o, such as a model includ- 
ing background emissivity only. The large fiux and intensity 
uncertainties account for the uncertainty of the beta func- 
tion parameters (core radius, beta and central intensity) and 
background value. 



3.5 



(T„ relation 



Figure 8 shows the location in the L — a„ plane of RzCS 052 
cluster, with all the z > 0.8 clusters for which we found in 
literature T,Lx and cr„ (Table 2). We ignored a few tenta- 
tive av determinations based on small number of velocities, 
because affected by large errors and by Eddington bias (de- 
tailed in appendix B and in sec 3.6). Literature values of Lx 
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Figure 8. Lx — o"u relation for literature clusters (open 
points) and for RzCS 052 (close point). The (red) triangle is 
CL1604+4304, i.e. the cluster originally claimed to be X-ray dark, 
at the revised value of the cry determination (see text for details). 



are converted to the [1-4] keV band rest-frame to minimize 
systematics. 

RzCS 052 (closed circle) turns out to be one of the 
faintest and least massive (lower velocity dispersion) clus- 
ters known at z > 0.8. In spite of the cluster being optically 
selected, RzCS 052 has an X-ray luminosity appropriate for 
its mass (velocity dispersion), or at least, data are compati- 
ble with the trend seen for similar, but X-ray selected, clus- 
ters. Therefore, although the X-ray luminosity of RzCS 052 
is modest, it is consistent with measurements for other high 
redshift clusters. This is reassuring for ongoing X-ray clus- 
ter surveys, which usually assume a Gaussian model for the 
scatter between mass and X-ray luminosity: RzCS 052 is 
not an example of a new class of clusters, massive but dark 
in X-ray, easily escaping the detection in X-ray surveys be- 
cause faint for their mass. Instead, RzCS 052 has an X-ray 
luminosity appropriate for its mass (velocity dispersion) and 
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Figure 9. Lx — Ov relation for literature clusters (close black 
points), as published by Fang et al. (2007), and DEEP2 groups 
and clusters, after our revision (red points and arrows). DEEP2 
groups and clusters are not longer underluminous for their veloc- 
ity dispersion. 

it is missed in the XMM Large Scale Structure because the 
low survey sensity at z ~ 1 for objects of RzCS 052 and 
obeying to the Lx — o"u relation. 

3.6 Are there known optically underluminous 
clusters or groups at high redshift? 

With the Bayesian tools described in the Appendix, we re- 
visit claims about the existence of underluminous X-ray 
clusters (i.e., clusters whose masses - from their velocity dis- 
persions - are too large for their X-ray luminosity - which 
is often an upper limit). We have seen (Appendix B) that 
in the case of C11604-I-4304 a Bayesian estimate of its veloc- 
ity dispersion, as well as the revised value of o\, published 
by Gal & Lubin (2004), make this cluster no longer an out- 
lier (i.e., underluminous) in the Lx — o^v relation. We show 
here that Bayesian estimates of velocity dispersions and X- 
ray fluxes cast serious doubts on the existence of the X-ray 
underluminous groups or clusters claimed in literature. 

Fang et al. (2007) study seven DEEP2 groups at 0.75 < 
z < 1.03. They derive an upper limit to the X-ray flux by 
considering only photons within an aperture of radius 30", 
a metric radius of 250 kpc which is considered to be typical 
for groups and clusters at high redshift. However, assuming 
j3 — 2/3 (also a typical value), the X-ray flux outside of their 
aperture, integrated to infinity, is 2.4 times larger then the 
fiux inside their aperture. The upper limit quoted by Fang 
et al. (2007) is therefore too small by a factor of 3.4. 

Fang et al. (2007) derive their velocity dispersion from 
3 to 6 galaxies. Because of Eddington (1940) biases, this is 
biased high: even symmetric errors move more low veloc- 
ity systems to high velocity dispersion than otherwise. The 
Bayes theorem allows to correct for the bias, being the Ed- 
dington correction built in the Bayes theorem (Appendix 
B). 

If these sources of error and biases are accounted for, 
the X-ray fluxes and velocity dispersions for the groups stud- 
ied by Fang et al. (2007) are perfectly consistent with the 
local Lx — cTv relation (Figure 9) . We note here that a sim- 
ilar argument can be made for the 'underluminous' CNOC 
groups claimed by Spiegel, Paerels, & Scharf (2007), based 



on velocity dispersions computed on just 3 or 4 velocities, 
and upper limits on the X-ray luminosity. 

Popesso et al. (2007) also claim that there exist X-ray 
underluminous clusters in the local universe, but their defl- 
nition of 'underluminous' depends on the data depth. Most 
of their underluminous clusters have normal X-ray luminos- 
ity for their mass, because they obey to the Lx — cr^ relation 
(see their figure 2d), and are called 'underluminous' because 
their are faint in their X-ray imaging. Deeper data would 
have classified them as 'normal'. The few remaining clus- 
ters are found to have a negative (unphysical) X-ray fiux 
and seem underluminous in the Lx — relation because 
they are plotted at an arbitrary value of Lx rather than as 
an upper limit. Such objects are not underluminous in the 
generally understood sense. 

To summarize, to our best knowledge there is no ev- 
idence for not even a single example of cluster (or group) 
of galaxies X-ray dim for its velocity dispersion (once all 
sources of errors are accounted for), all previous claims 
proven to be based on uncertain grounds. If there are un- 
derluminous clusters, they have not yet been convincingly 
discovered. 



4 HALO OCCUPATION NUMBER 

We now derive the HON for RzCS 052. This basically re- 
quires an estimate of the N galaxies within a specified lumi- 
nosity range and within a given cluster-centric radius and a 
mass estimate within the same spherical volume. 

The virial radius, r2oo, and mass, M200 are derived 
from the virial theorem from the measured cr„: r2oo = 1-04 
Mpc, and Khm = 4.0 1Q^*Mq. Adopting a Navarro, Frenk 
& White (1997) profile with concentration Cdm = 5 we 
derived rsoo = 0.69 Mpc and cluster mass within rsoo: 
M500 =2.9 10'" Mq. 

The number of galaxies is computed by integrating the 
luminosity function (LF, hereafter) within rsoo down to 
M* -\- 3. The latter is derived from the deeper VLT z' data 
presented in Andreon et al. (2007). Since the latter work 
count red galaxies only in their LF, we correct for the blue 
galaxies fraction (adopting the blue fraction measured in 
Andreon et al. 2007). We also correct for the distribution of 
cluster members in a 'cylinder' outside the cluster sphere, as- 
suming a NFW distribution with concentration Cgai = 3 (as 
in Lin et al. 2004, 2006). We opted for a Bayesian approach, 
because it simplify the computation of the uncertainty on 
the cluster richness fully accounting for uncertainties and 
covariances (neglected in past works) for all (Schechter and 
background) parameters. We found: A'500 = 41 ± 11 galax- 
ies. Fig. 10 compares RzCS 052 to the local richness-mass 
scaling, showing that it is within the local relation, although 
near the bottom-end of the distribution. 

Following Lin et al. (2006) we parametrize the evolution 
of N{m, z) as: 

N{M, z) = iVo (1 + zf {M/MoY (1) 

where No — 56 and Mo = 2 lO'** are normalization fac- 
tors of the relation, and s = 0.84 is the slope of the local 
relation derived by Lin et al. (2004) for their local cluster 
sample. We can rewrite this as below, to emphasize the evo- 
lutionary terms: 



8 Andreon et al. 




log [M] redshift 

Figure 10. Left panel: Number of galaxies as a function of cluster mass. Open points mark local (z < 0.1) clusters, the closed (red) 
point mark RzCS 052. The line is the Lin et al. (2006) fit to local clusters. Right panel: Evolution of the N — M relation. Points are as 
left panel. The shaded (yellow) region mark the 68 % confidence region derived by Lin et al. (2006) at 2: < 0.9. 



(1 + zy° = 



No{M, 



No (M/Mo) 



(2) 



where we have added the subscript to emphasize that 
we are now talking about the observed (or maximum like- 
lihood) values. In order to estimate 7 we can just look at 
the dependence of the r.h.s of eq 2 with (1 + 2), as shown 
in the right-hand panel of Fig. 10. The value observed for 
RzCS 052 is a bit farther away in redshift than the range 
probed by Lin et al. (2006), and is outside their 68 % in- 
terval on 7, shown as shaded (yellow) area. A zero value 
would implies that the way galaxies populate cluster-scale 
halos at z = has not changed from z = 1. The size of 
the error bar on RzCS 052 is comparable to the error on 7 
(the width of shaded region in figure at z = 1.0), indicat- 
ing that RzCS 052 alone carries comparable information to 
all the high redshift clusters studied by Lin et al. (2006). 
Therefore, the data for RzCS 052 suggest a mild evolution, 
with the caveat that the scatter around the mean relation 
is large (see left-hand panel) and our result should be taken 
as tentative. We note, however, that our measurement is 
more direct than Lin et al. (2006): we include an estimate of 
the characteristic luminosity M* and faint-end slope a from 
our data, whereas Lin et al. (2006) assumed them for lack of 
data, and we measure mass and reference radius, rsoo, from 
the virial theorem without assuming that they scale with 
X-ray temperature and evolve self-similarly, as assumed by 
Lin et al. (2006) for lack of direct measurements. The mild 
difference seen in the right panel of Fig 10 may indicate a 
possible break in the (assumed) self-similar evolution of the 
scaling between temperature and radius or mass at 2 ~ 1. 

The parameter 7 is the (logarithm) derivative of the 
redshift dependence of the number of galaxies per unit clus- 
ter mass, i.e. of the galaxy merging rate in appropriate units. 
Lin et al. (2006) and our results agree that the number of 
galaxies per unit cluster mass has increased (us) or stayed 
constant (Lin et al. 2006) since z = 1. Therefore, both stud- 
ies directly show that no intense merging activity of galaxies 
has been ongoing in clusters in the last 7 Gyr. 



5 SUMMARY 

We have identified a distant cluster from a modified red- 
sequence method and followed it up spectroscopically. RzCS 
052 is a richness class 3 cluster at z = 1.016 with a velocity 
dispersion of 710 ± 150 km s"^ and an X -ray luminosity of 
0.68 ±0.47 X 10^* ergs s~\ 

In spite of its optical detection, RzCS 052 obeys to the 
high redshift Lx ~ o"u relationship as other X-ray selected 
clusters to the high redshift Lx — o"„ relationship, whereas 
in principle variations in the dynamical state of the clusters 
or in the thermal history of the intracluster medium may 
have moved it away from the Lx — o"u relation. 

Analysis of the N — M scaling shows that RzCS 052 
has the right number of galaxies (actually, a bit less) that it 
should have for its mass, ruling out intense merging (among 
galaxies) activities in clusters from 2 = 1 to today, in agree- 
ment with Lin et al. (2006). 

We present a Bayesian approach to measuring cluster 
velocity dispersions (most useful for sparsely sampled data 
and in presence of a background) and X-ray luminosities or 
upper limits (essential in the case of poorly determined pa- 
rameters). Critical re-analysis of the data of clusters/groups 
claimed to be outliers of the Lx — o"u relationship leads to 
conclude that there are no known thus far examples of clus- 
ters X-ray underluminous for their velocity dispersion. The 
above result is quite reassuring for ongoing X-ray surveys: 
there is thus far no example of cluster missed because an 
anomalous Lx for the cluster mass. 
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Table 2. Luminosity and velocity dispersion of clusters at 2 > 0.8 



Name 


z 


log Lx [1-4] 


ref 




N 


ref 


RXJ1716+6708 


0.813 


44.68 ±0.03 


12 


1522 ± 180 


37 


1 


RXJ1821. 6+6827 


0.816 


44.62 ± 0.01 


2 


775 ± 122 


18 


2 


MS1054-30321 


0.830 


44.91 ±0.04 


12 


1153 ± 80 




3 


RXJ0152-1357S 


0.830 


44.42 ± 0.02 


12 


737 ± 126 


18 


4 


RXJ0152-1357N 


0.835 


44.59 ±0.03 


12 


919 ± 168 


16 


4 


RzCS 530 


0.839 


44.20 ±0.09 


5 


780 ± 126 


17 


5 


1WGA1226+3333 


0.890 


45.17 ±0.01 


12 


997 ± 245 


12 


6 


C11604+4304 


0.900 


44.00 ±0.06 


14,7 


962 ± 141 


67 


7 


RzCS 052 


1.016 


43.83 ±0.37 


this work 


710 ± 150 


21 


this work 


RXJ0910+5422 


1.106 


44.00 ± 0.05 


12 


675 ± 190 


25 


8 


RXJ1252-2927 


1.237 


44.37 ± 0.07 


12 


747 ± 79 


38 


9 


LynxW 


1.270 


43.70 ±0.23 


12 


650 ± 170 


9 


10 


1WGAJ2235.3 


1.393 


44.57 ± 0.05 


11 


762 ± 265 


12 


11 



RzCS 530 is also known as XLSSC 003; References: 1: Gioia et al. (1999); 2: Gioia et al. (2004a); 3: Gioia ot al. (2004b); 4: Demarco et 
al. (2005); 5: Valtchanov ct al. (2004); 6: Maughan et al. (2004); 7: Gal & Lubin (2004); 8: Mei et al. (2006); 9: Demarco et al. (2007); 
10: Stanford et al. (2001); 11: MuUis et al. (2005); 12: Ettori et al. (2004); 14: Lubin et al. (2004). 
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APPENDIX A: MIXTURE MODELLING OF 
INHOMOGENEOUS PROCESSES FOR THE Lx 
AND RICHNESS ESTIMATES 

We want to rnoasuro a structured (i.e. not constant in space) 
Poisson signal in presence of a background. Wo assume that 
the background (photons, galeixies, etc.) distribution is an 
homogeneous (i.e. the intensity is independent of position) 
Poissonian random process, whereas the cluster contribu- 
tion is an inhomogenous Poissonian random process whose 
intensity is given by an 7(r) radial profile. Provided that 
quantities are Poisson distributed, it does not matter if we 
are talking about X-ray photons (as in sec 3.3), or galaxies 
(as in sec 3.4) or something else. 

Let us call 6 the (unknown) set of parameters of the 
function /(r). Simple algebra shows that the likelihood func- 
tion £(61) =p{ri\I{r)) is 

m = Jlu;{n)I{n) e-X,'^^'"^'^'"^ (Al) 

i 

where is the solid angle. The expression can be sim- 
plified somewhat by noting that the infinitesimal solid angle 
at Ti, uj{ri), is independent of the parameter d and therefore 
can be dropped. There are no limitation on the complexity 
of the shape of the solid angle. O also encodes the relative 
efficiency of the different parts of the instruments (for ex- 
ample, the different efficiency of off- and on- axis response). 

Combined with prior probability distributions for the 



parameters, this likelihood function yields, via the Bayes 
theorem, the posterior distribution for the function pa- 
rameters d, given the data. Marchov-Chains Monte Carlo 
(Metropolis 1953) with a Metrojjolis (1953) sampler is used 
to sample the posterior. The chain provides a sampling of 
the posterior that directly gives credible intervals for what- 
ever quantity, either for the parameters 9 or any derived 
quantity such as total richness (or flux): for am interval at 
the desired credible level it is simply matter of taking the 
interval that includes the relevant percentage of the sam- 
plings. Credible intervals (yellow areas in Fig 4 and Bl) are 
computed in that way. Upper limits may be determined in 
the same way as fluxes for detection, i.e. by specifying the 
credible interval we arc interested in. 

The function 7(r) can be whatever function positively 
defined and having a finite integral. In this paper we use a 
modified (3 function: 

I{r) oc [(1 + r/rc)']-'^+'/' + hkg (A2) 

were we have accounted for a constant background, hkg. 
By choosing a more complex background function, as in An- 
dreon (2006b), we obtain the aimed mixture modeling of two 
inhomogenous Poissonian processes. 

Al Fallacies of the usual measurements of Lx 
upper limits 

While our way of determining fluxes, richness and their er- 
rors, as well as Lx upper limits is unusual in our astronom- 
ical context (but the standard approach in other fields of 
astronomy and in statistics), wo were obliged to introduced 
it because previous approaches arc unsatisfactory when an 
important parameter has a large error or is undetermined. 

A common assumption of many determinations of upper 
limits to the X-ray fiux from a cluster is that the object flux 
is fully inside a given aperture or the object core radius and 
(3 are known. However, this is a dangerous assumption: if the 
object is undetected, its extent, core radius and (3 axe not 
constrained. If, for example, the object is much larger than 
assumed or /3 is small (and data tell nothing about that, 
being the object undetected) the assumption has important 
consequences. For example, it is sufficient to assume that the 
'underluminous' groups of Fang et al. (2007) have 'typical' 
core radii to make their Lx compatible with the local Lx — 
av relation (and the groups no longer underluminous). 

A first step in the right direction is to correct for the 
fiux outside the aperture, but this assume to known the 
unknown: when rc and (3 are unmeasured, or are very poorly 
determined, we cannot assume them as perfectly known and 
we cannot make inferences dealing with quantities strongly 
depending on the poorly determined parameters, such the 
location of clusters in the Lx — cTv relation. The scientific 
method does not suggest to hope to have taken, by good 
chance, the correct value of an unknown parameter (as rc) 
when it strongly affects the result. 

However, this is exactly the kind of problem where 
Bayesian approaches are most valuable. Bayes' theorem al- 
lows us to infer the value of a quantity (in this case Lx) 
in the presence of a nuisance parameter (core radius or (3) 
whose value is unknown but whose value affects the measure- 
ment of the quantity. Assuming a single value for nuisance 
parameters artificially collapses the error ellipse along one 
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-2000 2000 
rest — frame velocity 

Figure Bl. Velocity histogram for a simulated dataset com- 
posed of 500 member galaxies drawn from a Gaussian having 
true av = 1000 km/s and 500 interlopers drawn from an uniform 
distribution. The red curve is the cluster velocity estimate de- 
rived from robust statistics {a^ = 1400 km/s), whereas the blue 
curve and the shaded yellow region show the Bayesian estimate: 
av = 940 it 85 km/s. See text for details. 



(or more) axes and leads to an incorrectly small error bar 
and to call outlier something that instead is fully compatible 
with the model. The sum rule of probability prescribes to 
marginalize (average over) nuisance parameters, not to keep 
them fixed. 

Other authors determine upper limit to the X-ray flux 
mistakenly taking the maximum likelihood estimate of sam- 
pling theories, total — background, for the true value of the 
net flux. While total — background is allowed to be negative, 
the true value of the net flux cannot be. These two quan- 
tities differ when the net flux is comparable to background 
fluctuations (e.g. appendix B of Andreon et al. 2006). 

Another common way to compute the upper limit of the 
X-ray flux is by measuring the fluctuations of background 
counts. While this number is interesting in its own right 
and has the appealing property that it becomes smaller and 
smaller with lower and lower background fluctuations, it is 
measuring something different than the X-ray flux. In fact 
this quantity is a p— value, i.e. a measure of how frequently 
one observes larger background fluctuations under the null 
hypothesis that no (cluster) signal is there, which differs 
from how probable a signal can be there without detecting 
it. A pedagogical astronomer oriented explication of the dif- 
ference of the two concepts is presented in Andreon (2008) . 



APPENDIX B: MIXTURE MODELING FOR av 
MEASUREMENTS 

We want to measure the scale (dispersion) of a distribution 
(say, of velocities) , knowing that the sample is contaminated 
by the presence of interlopers, but without the knowledge of 
which object is an interloper. The main idea is not identify 
or de-weight interlopers in the scale estimate, but account 
for them statistically, precisely as astronomers do with pho- 



tons when estimating the flux of a source in presence of 
a background. The small size of astronomer samples (e.g. 
of cluster galaxies with known velocity) makes the asymp- 
totic properties of frequentist estimators never reached in 
real life experiments and oblige us to look for a solution in 
the Bayesian paradigm. 

Here, we assume that data come from two populations: 
background galaxies, whose distribution is assumed to be an 
homogeneous (i.e. the intensity is independent on v) Poisso- 
nian random process, and cluster galaxies, whose distribu- 
tion is assumed to be a Poissonian process whose intensity is 
given by a Gaussian. The likelihood is given by eq Al, with 
changes of variable names: Q continues to be fi, but it is eas- 
ier to understand it if we call it Av, the (velocity) range over 
which velocities are considered (say, ±5000 km/s from the 
cluster preliminary velocity center). Q is more appropriate 
than Av as it accounts for intervals of complicated shape; r 
in Appendix A is now v, and I{v) is given by the sum of a 
Gaussian and a constant, with unknown weights, Ndus and 
Nbkg (respectively). Finally, each measured velocity v has 
an uncertainty a^- Therefore, I{v) reads: 



I{v) 



+ 



N, 



bkg 



Av 



(Bl) 



Most literature estimates of cluster velocity dispersions 
are based on the family of estimators presented by Beers, 
Flynn & Gebhardt (1991). However, in the presence of a 
background and interlopers, with sparsely sampled data, a 
Bayesian estimator may be more appropriate. 

A 'real life' example may suffice. The cluster 
C11604-f 4304 was regarded as unusually X-ray dim for the 
large mass estimated from a sample of 27 redshifts o\, = 
1226lig km s~^ (Postman, Lubin & Oke 2001). However, 
the Bayesian method returns (t„ = 1022 ± 570 km 
(posterior mean and standard deviation), which no longer 
makes the cluster X-ray underluminous and has a more re- 
alistic error bar. A larger sample of redshifts for this clus- 
ter, from Gal & Lubin (2004) revises the original estimate 
to Ov — 962 ± 141 km s~^, in good agreement with the 
Bayesian estimate. With this value, the cluster is no longer 
X-ray underluminous. 

Let us now consider a simulated 'cluster' composed of 
500 galaxies distributed in a Gaussian with cr„ = 1000 km 
s~^ and superposed over a background of 500 uniformly dis- 
tributed (in velocity) interlopers. The large sample size has 
been adopted to leave data to speak by themselves. Apply- 
ing the methods of Beers et al. (1991) yields S\, — 1400 
km s~^ which is an excessively large estimate of a„ (and 
hence of mass), as also visible in Fig. Bl by simple inspection 
(compare the red curve and the histogram). The Bayesian 
posterior mean is = 940 ± 85 km s"'^ which is closer to 
the 'true' value (blue curve with shading). This simulation 
shows that the amplitude of bias of the Beers et al. estimator 
is systematic (i.e. it is present even for a large sample), and 
it is actually independent on the sample size, provided the 
relative fraction of cluster and interlopers is kept, although 
harder and harder to note as the sample size decreases be- 
cause the estimator variance increases and dominates the 
scatter. 

We now assess the sensitivity to model assumptions. 
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0,5 




Av [1000 km/s] 

Figure B2. Perturbed velocity distribution (solid line), given by 
eq. B2, and a Gaussian with identical first two moments (dashed 
blue line). The former is used to generate hypothetic data, the 
latter is assumed to estimate cr„. 



the data does not tightly constraint the aimed quantity) , be- 
cause, as well known to astronomers, if there are many more 
low velocity systems than high velocity dispersion systems 
the observed value (i.e. the maximum likelihood value) is a 
biased estimate of the 'true' value. As point out by Jeffreys 
(1938), the Bayes theorem quantifies the bias, and we used 
it for computing the correction to Fang et al. (2007) veloc- 
ity dispersions. Specifically, we assume a logarithmic slope 
of —0.6 for the prior and we follow appendix A of Andreon 
et al. (2006), because Fang et al. (2007) do not publish in- 
dividual velocities for their systems. 

Beers, Flynn & Gebhardt (1991) scale estimators work 
correctly in many cases, as shown in their paper. In these 
cases, the Bayesian approach returns similar numbers. We 
have shown, however, that in frontier-line cases, i.e. in the 
presence of an important background, or with sparsely sam- 
pled data, the Bayesian method returns better behaved 
quantities, less biased and less noisy. 



Lets suppose that cluster substructure perturbs the velocity 
distribution, that we now assume to be described by 

p(«) cxe''/i°»''(l + e2-^5"/™°)-^ (B2) 

depicted in Fig B2 (solid line). The function has first 
and second moments (mean and dispersion) equal to —460 
and 1130 km s~^, respectively. We simulate 1000 (virtual) 
clusters of 25 members each (and no interlopers) drawn from 
the distribution above (eq B2), but we compute the velocity 
dispersion using eq. Bl, i.e. with a likelihood function ap- 
propriate for members drawn from a Gaussian. The mean of 
found posterior means is = 1140 km s~^ (vs. the 1130 km 
s~^ input value) with a standard deviation of 185 km s~^. 
The mean error uncertainty (posterior standard deviation) 
is 163 km s~^, close (as it should be) to the scatter of the 
posterior means. The uncertainty has a negligible scatter, 
18 km s~^, indicating the low noise level of each individual 
uncertainty determination, four time lower than the scatter 
of the uncertainty of the biweight estimator of scale (70 km 
s~^), that instead shows values as small as 73 km s~^ and 
as large as 865 km s~^ for data that are supposed to give a 
unique, fixed, value of uncertainty. 

As more difficult situation, we now consider a sample 
drawn, as before, from a distribution different from the one 
used for the analysis, but furthermore ~ 50% contaminated 
by interlopers and consisting of half as many members: 13 
galaxies are drawn from the distribution above (eq. B2), 
superposed to a background of 12 galaxies, uniformly drawn 
from ±5000 km s~^. The mean of found posterior means is 
cr„ = 1160 km s~^ (vs. the 1130 km s~^ input value). The 
mean error uncertainty is 390 km s~^, with a low (80 km 
s~^) scatter. The biweight estimator returns, on average, a 
strongly biased estimate o\, — 2135 km s^^. 

The Bayesian determination of the cluster velocity dis- 
persion already embodies the correction for the Eddington 
bicis: the prior (i.e. the number distribution of objects hav- 
ing does matter when the likelihood is shallow (i.e. when 



